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Abstract 



We show that every n-vertex cubic graph with girth at least g have domination 
number at most 0. 299871™ + O (n/g) < 3n/10 + O (n/g). 

u 

ji ■ 1 Introduction 

-i— > 

The notion of a dominating set is a classical notion in graph theory with a large amount 
of literature associated with it. For the sake of completeness, let us recall that a set D of 
vertices of a graph G is dominating if every vertex of G is contained in D or has a neighbor 
in D and the domination number 7(G) of G is the smallest size of a dominating set of G. 
\q \ In this paper, we study the domination number of cubic graphs, i.e., graphs where every 
vertex has degree three. 

In 1996, Reed [9] conjectured that every n-vertex connected cubic graph has domination 
number at most [Vi,/3~|. Though the conjecture turn out to be false [HE], the conjecture 
becomes true with additional assumption that the cubic graph has girth at least g, i.e., it 

O 



has no cycles of length less than g. The first results in this direction are the bounds on 
the domination number of ^| + j n of Kawarabayashi et al. [2 J for bridgeless n- vertex 

cubic graphs with girth at least g for g divisible by three and ^| + n of Kostochka and 

Stodolsky [5] for all n- vertex cubic graphs with girth at least g. The magic threshold of n/3 
was first beaten for cubic graph with large girth by Lowenstein and Rautenbach [6] who 
showed that every n-vertex cubic graph with girth at least g > 5 contains a dominating 



*This research was supported by the Czech-Slovenian bilateral project MEB 090805 (on the Czech side) 
and BI-CZ/08-09-005 (on the Slovenian side). 

^Institute for Theoretical Computer Science (ITI), Faculty of Mathematics and Physics, Charles Uni- 
versity, Malostranske namesti' 25, 118 00 Prague 1, Czech Republic. E-mail: kral@kam.mff.cuni.cz. 
Institute for Theoretical computer science is supported as project 1M0545 by Czech Ministry of Educa- 
tion. This research was also supported by the grant GACR 201/09/0197. 

■^Department of Applied Mathematics, Faculty of Mathematics and Physics, Charles University, Mal- 
ostranske namesti 25, 118 00 Prague 1, Czech Republic. E-mail: peskoj@kam.mff.cuni.cz. 

^Department of Applied Mathematics, Faculty of Mathematics and Physics, Charles University, Mal- 
ostranske namesti 25, 118 00 Prague 1, Czech Republic. E-mail: janv@kam.mff.cuni.cz. 



1 



set of size at most ^ ^ + iff^J n ~ 0.325926n + O {n/g) The bound was further improved 

by Rautenbach and Reed [8] to 0.321216n + O (n/g). 

We further improve these bounds and manage to lower them below the 3n/10 threshold. 
Our main result is the following: 

Theorem 1. Let G be an n-vertex cubic graph with girth at least g. The domination 
number of G is at most 0.299871n + O {n/g) < 3n/10 + O (n/g). 

At this point, we remark that numerical computations involved in our argument were done 
using a computer (but the rules given in Figure [3] were generated by hand) though the 
whole proof can be easily verified to be correct without computer assistance. 

Before we start the exposition of our proof, let us mention a connection to random cubic 
graphs. A random cubic graph almost surely contains only a bounded number of cycles of 
length less than g for every fixed integer g and thus the value of the domination number of 
a random cubic graph should indicate how tight our result could be. It is known that the 
domination number 7 of a random cubic graph is almost surely at least 0.2636n [7] and at 
most 0.2794n [1] which should indicate space for further improvements. 



2 Proof 

We first provide a general overview of our method and illustrate it on a small example. At 
the end, we then apply the method in a setting yielding Theorem HJ 



2.1 Overview 

In this subsection, we explain main ideas of our method and we later provide necessary 
technical details related to it. Fix a cubic bridgeless graph G with girth at least g and also 
fix an integer K which determines the number of levels as defined later. Consider a 2-factor 
of G (which exists by the Petersen theorem) and decompose each cycle of the 2-factor into 
vertex-disjoint paths Pi, . . . ,P# with the number of vertices between g/^K and g/2K (this 
is possible since the length of each cycle of the 2-factor is at least g) . The vertices of the 
paths Pi, . . . , Pi are considered to be ordered from one end of the path towards the other; 
a mate of a vertex v is the neighbor of G not adjacent to v on the cycle of the 2-factor. 

A dominating set D of G will be given by a labeling of vertices of G we construct. 
Each vertex of G will be assigned an input label which is one of the symbols +, x , • and 
o, and an output label which is one of the symbols ©, Cg> and 0. The dominating set D 
will contain the vertices with input label + or output label <S> (as well as several others, see 
the next subsection for an exhaustive definition). We explain the intuitive meaning of the 
labels later. 

Split the paths P^ . . . ,P k into K sets Vi, ... , Vk including each path to a single set 
randomly uniformly and independently of the other paths. The sets Vx, . . . , Vk are referred 
to as levels and vertices on paths in Vi are said to be on the level i. First, the vertices 
contained in paths of V\ are assigned labels, then those in paths of V2, etc. Let P be 
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Figure 1: An example of a (correct) set of rules. 

a path included in Vi and assume that the vertices of the paths in V\ U • • - U T-^-i have 
already been labelled. The input label of a vertex of P is: 

— the symbol + if its mate is on a path in V\ U ■ • ■ U V%-\ and its output label is 0, 

— the symbol x if its mate is on a path in V\ U • ■ ■ U Vi-i and its output label is 0, 

— the symbol • if its mate is on a path in V\ U • • • U V%-\ and its output label is 0, 

— the symbol • if its mate is on a path in Vi, and 

— the symbol o if its mate is on a path in 7-Vn U • • • U TV 

The output labels are assigned in blocks using rules. Each rule is a pair of a sequence of 
input symbols and the symbol ? which represents a wild-card and a sequence of output 
symbols. The lengths of the two sequences will always be the same. We usually write an 
arrow between the two sequences, e.g., one of the rules can be +? — > 00. Naturally, the 
sequence of input symbols and ? is called the left-hand side of the rule and the sequence of 
output symbols the right-hand side. Finally, if a — > r is a rule, we use Cj is the i-th symbol 
of a and 7$ is the i-th symbol of r. An example of a set of rules is given in Figure [TJ 

We look for a rule whose left-hand side matches the input labels of vertices at the 
beginning of P. All the considered rule sets will have the property that such a rule is 
unique. The right-hand side of the rule then determine the output labels. We then move 
after the vertices with output labels assigned and again look for a rule whose left-hand side 
matches the input labels of vertices without output labels. After finding a suitable rule, 
the next group of vertices is assigned output labels, we move after them and continue in 
this way until we reach the end of the path P. Note that it can happen that we are left 
with few vertices without output labels at the end of the path P — this will be handled in 
the next subsection. 

The meaning of output labels is the following: if a vertex v is labelled with <g>, then v 
should be included to D. If a vertex v is labelled with ©, then its mate should be included 
to D (this output label can only be assigned to vertices with input labels o, i.e., with 
mates on higher levels). Finally, if v is labelled with 0, then it is dominated by one of its 
neighbors on its path or by its mate on a lower level. The input labels can be interpreted 
as follows: the label + represents that v is included to D, the label x represents that its 
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mate on a lower label is included, the label • represents that neither v nor its mate has 
yet not been included to D and the mate is on a lower or the same level and the label o 
represents that a mate is on a higher level. 

A set of rules is called correct if for every rule a —>■ t the following holds: 

— if Oi is • or o and neither nor o~ i+1 (if they exist) is +, then n is © or one of the 
symbols Tj_i, Tj and r i+ i is ®, and 

— if ^ is ©, then <7j is o. 

Intuitively, if a set of rules is correct, then the set D containing the vertices with input 
label + or output label ® is always dominating and the label © can be only assigned to 
vertices with mates on higher levels. The set of rules given in Figured] is correct. 

2.2 Analysis and adjustments 

In this subsection, we provide further details on the labeling procedure and analyze it. The 
first thing to cope with is the fact that a cubic graph need not to have a 2-factor. This is 
handled in a way analogous to that used in |8j. For a collection of vertex-disjoint paths, a 
vertex is covered if it is contained in one of the paths. 

Lemma 2. For every K > 2, g and n-vertex cubic graph G with girth at least g, there 
exists a collection of at most 3+ 2 g K n vertex disjoint paths in G with the number of vertices 
less than g/2K that covers at least n — O (n/g) vertices of G. Moreover, all vertices on the 
paths can be grouped into pairs of vertices adjacent through an edge not contained in the 
paths. 

Proof. We only sketch the proof as it essentially follows the lines of the argument in [8]. 
Let S be a set of vertices of G. The number of odd components of G \ S with a single edge 
to S is at most n/g (since each such component contains a cycle; in fact their number can 
be bounded by n/2 e<y9 \ but we do not need this finer estimate here). Hence, the Tutte- 
Berge formula implies that G has a matching M covering all but at most n/g vertices of 
G. Remove the vertices not covered by M and the edges of M. In this way, we obtain a 
collection of at most 3n/2g paths and cycles, each of them with length at least g. 

We split the cycles into paths with the number of vertices at least g/4K and less than 
g/2K. This is clearly possible. Similarly, the paths with at least g/2K vertices are split. 
Paths with less than g/2K are preserved and there is at most 2>n/2g such paths. Altogether, 
the final collection of paths consists of at most AKn/g new paths (as each of them has at 
least g/^K vertices) and at most 2>n/2g paths of the original collection. □ 

Fix a cubic graph G and a collection of vertex disjoint paths Pi, . . . , as described 
in Lemma [2l We say that a vertex w is 1-close to another vertex v if w lies on the same 
path as the mate of v (in particular, w can be the mate of v). Note that the relation of 
being "1-close" is not symmetric. For i > 1, we say that w is i-close to v if w is 1-close to 
v or w is (i — 1) -close to a vertex w' where w' is not the mate of v but lies on the same 
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path as the mate of v. Observe that if w is i-close to v, then the distance between v and 
w is less than gi/2K. Because of this, there is no vertex w that is .fT-close to two different 
vertices v and v' lying on the same path (if such w existed, then the paths from v and v' 
to w witnessing that w is K-close and the subpath between v and v' would form a cycle of 
length less than g) . The notion of being K-close will play a crucial role in our argument 
that certain events are independent later in the proof. 

For each path P among P 1 , . . . , P^, choose randomly uniformly and independently of 
other paths an integer i, 1 < % < K, and include P in the set TV We now recursively define 
the probabilities Pi{+), Pi(x), Pi(*) and Pi(°) which represent probabilities that a vertex 
of a path in Vi has a certain input label, the probabilities (&(©), <7i(<8>) and g«(0) that it 
has a certain output label and g°(©), <7°(<S>) and ?f(0) that it has a certain output label 
conditioned on input label being equal to o. These numbers represent the probabilities in 
the "ideal" case, which includes the assumption that all the labels are "independent" of 
each other and all the paths in Vi are infinite, and we will have to argue that they can also 
be applied to our labeling procedure. 

Let us start with estimating the input probabilities. A mate of a vertex v of a path in 
Vi is on a path of the lower level or the same level with probability i/K, thus 

Pi(o) = 1 - ± ■ (1) 

Since the mate is on a path of any fixed level with the same probability equal to 1/K, we 
also obtain 

Pi(+) = E^?» (2) 

j'=i 

Pi(x) = E^?» (3) 
j'=i 

= ^ + £^ (0) (4) 

Note that if « = 1, then Pi(+) = 0, p^(x) = and = 1/K. 

Once we have determined the probabilities of input labels, we can compute the proba- 
bilities of output labels (assuming that all input labels on a single path are independent of 
each other). If 1Z is a set of rules and a — > r is one of the rules, we use |<r| for the length 
of a, \t\ for the length of r (which are the same in each individual rule) and ar(x,y) for 
the number of symbols y in r on positions of symbols x in a (with x = ? or y = ? being 
the wild-card), e.g., if o = o« + and r = ® 0, then |er| = \r\ = 3, crr(?,©) = 2 and 
ot(o,o) = 1. If a" is a sequence of input symbols, then Pi(cr) denotes the probality that 
the sequence a appears on |er| consecutive vertices as their input symbols, i.e., 

M 

pm) = n^^') • 
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Finally, \\7Z\\ is the maximum length of a rule in TZ and TZ is the expansion of 7Z, the 
set of rules obtained from TZ by replacing each rule containing the wild-card symbol(s) ? 
with rules for all possible choices of the value of ?, e.g., a rule containing two symbols ? is 
replaced with 16 rules. 

The probabilities of output labels of vertices are then given by the following formulas: 



Qi 



= ( 6 ) 
?*(©) = £_^rW,©) (7) 

where 

Let us derive (J5D — (ED formally. Consider a path of length L on level i. Let be the 
expected number of rules of TZ matched and E a the expected number of times the rule 
cr — > t is matched. In the analysis that follows, it is convenient to think of the considered 
path as a sufficiently long part of an infinite path which is consistent with our arguments 
presented later. When we start matching, the probability that the rule a — > r is matched 
is equal to Pi(cr). Hence, the expected number of times E a the rule o — > r is matched is 
Pi(a)E. By the linearity of expectations, we also have that 



£ EM = L. 



By the definition of Qi, we obtain that E = L/Qi and E a = Lpi(a)/Qi. Summing over 
all rules a — > r G TZ, we obtain the expected numbers of output labels of the vertices of a 
path, e.g., the expected number of the output label © is 

The obtained quantities after dividing L represent the corresponding probabilities as given 
in (ED-©. 

Similarly, the probabilities conditioned on the appearance of the input symbol o are 
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Figure 2: The probabilities given by f|T])- f|T0|) for the set of rules from Figured] and K — 5. 
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The just defined quantities for the set of rules given in Figure [1] and K = 5 are given in 
Figure El 

In the ideal case, we understand each of the paths in Vi as an infinite path where the 
probability that any particular vertex is matched to the first symbol in one of the rules of 
1Z is the same. This probability denoted as rn, is equal to 



E_Pi(<r) 



1 



Y, _W\Pi( cr ) Qi ' 

This formula can be derived by considering a path of sufficiently large length L and dividing 
the expected number of rules matched (which we computed earlier) by L. Similarly, we 
define r^- to be the probability that a vertex is matched to the first symbol of a rule in 1Z 
conditioned that none of the j — 1 previous vertices has been matched to the first symbol 
of a rule. Clearly, r^j = for j > \\7Z\\ = ||72.||. As an example, observe that n i = 1/3, 
r i,2 = 1/2 and r\ 3 = 1 for the set of rules from Figure [TJ 

We now describe the actual procedure used to label the vertices of paths and give all 
details to cope with technical difficulties which we omit in Subsection 12.11 The paths are 
labelled from the first level, i.e., we start with paths in V\, continue with those in V2, etc. 
In each set Vi, the paths are labelled in arbitrary order. 
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Consider a path P from V%. The input label of a vertex v is determined as described 
in Subsection 12. 1[ Next choose randomly a number £ between 1 and ||7?.|| whose value 
is equal to £ with probability rj/]X-Ii(l ~~ r «,i)- Add auxiliary \\1Z\\ — £ vertices at the 
beginning of P and label each of them with x with probability Pi(x) and add auxiliary 
1 1 1Z 1 1 vertices at the end of P and label each of them with x with probability Pi(x). Finally, 
apply the labeling procedure as described in Subsection 12.11 (starting with the £ auxiliary 
vertices) and when finished, discard the £q + \\TZ\\ auxiliary vertices. In this way, all vertices 
of P are assigned output labels. 

Let us analyze the actual probabilities that a vertex in V% has a certain input and output 
labels. We claim that a vertex of a path P of V% is assigned an input label x G {+, x, •, o} 
with probability Pi(x) and an output label x G {©, 0, 0} with probability qi{x). Moreover, 
the probability of a vertex v getting a certain input label depends only on the labels and 
the levels of vertices that are i-close to v (and thus the input labels of vertices on P are 
mutually independent random variables). Because of the addition of ||7?.|| — £$ auxiliary 
vertices at the beginning of P, where £q was chosen as described earlier, each vertex of P 
has the same probability of being the first, second, etc., in a rule of 1Z applied to P. Hence, 
the probability that its output label is y G {©, 0, 0} is equal to qi(y). 

We now analyze probabilities of input and output labels. Let v be a vertex of a path 
P in V%. Since the mate v' of v is on a higher level with probability 1 — i/K, the input 
label of v is o with probability 1 —i/K. With probability 1/K, the mate v' is on the 
same level and the input label of v is • in this case. With probability (i — 1)/K, the mate 
v' is on a lower level and the input label of v is determined by the output label of v' in 
this case; the probability that the output label of v' is y G {©, 0, ©} is q°(y) where j is 
the level of v' . As the output label of v' depends only on the levels and labels of vertices 
(i — l)-close to other vertices on the path of v', the output label of v' depends only on the 
vertices i-close to v (and the level of its mate). Hence, the input label of v depends only 
on the labels and the levels of vertices z-close to v. In particular, input labels of all the 
vertices of P are mutually independent. Since the probability of v being the first, second, 
etc. in a particular rule of 1Z applied to it is the same because of padding with £ auxiliary 
vertices, the probability that the output label of v is y G {©, 0, 0} is equal to qi(y) and 
the probability that the output label is y conditioned by its input label being o is q°(y)- 

Consider a labeling of the vertices of G constructed in the just described way. The 
dominating set D for a graph G is formed by the following vertices: 

— vertices not covered by the paths in V\ U • • • U Vk, 

— vertices with input label +, 

— vertices with output label 0, and 

— the first and the last vertex of each path in V\ U • • • U Vk- 

Let us verify that D is a dominating set: a vertex v not covered by the paths in V\ U ■ • • UVk 
is in D. Vertices on paths V\ U • • • U Vk with input label different from + and output label 
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different from (g) are dominated either by their mates or by their neighbors on paths with 
auxiliary vertices (assuming the set of rules is correct). However, since the auxiliary vertices 
were discarded, the first and the last vertex of each path may not be dominated in this 
way — this has been repaired by adding them to D (regardless their labels). 

It remains to estimate the expected size of D. The expected size of D is equal to the 
sum of the expected number of vertices with output label <8> or © (note that each vertex 
with input label + has a mate with output label ©), the number of vertices not covered by 
paths (which is at most 0(n/g)) and twice the number of paths because of the inclusion 
of their first and last vertices to D (this number is also 0(n/g)). The probability that a 
vertex v has output label © is ag ^ e path containing v is included to any 

of the K levels with the same probability. Similarly, its output label is © with probability 

K 

We summarize the results presented in this section in the following lemma: 

Lemma 3. Let TZ be a correct set of rules, K > 2 an integer and gj(ffi) and qi(<2)) quantities 
determined using /fi])-/ f70j) . If the procedure described in this subsection given by the set of 
rules TZ and the integer K is applied to any cubic n-vertex graph with girth at least g, then 
it produces a dominating set with expected size equal to 



££i(ft(©)+ft(®)) n + 



7? 



K \g, 

Note that already the simple set of rules given in Figure [T] for K = 5 yields by Lemma [3] 
that the domination number of a cubic graph with girth at least g is at most 0.313972n + 
0(n/g), an improvement of the bound of Rautenbach and Reed from [8]. 

2.3 Finale 

We apply Lemma [3] for a suitable correct set TZ of rules. This set of 79 rules can be found 
in Figure [31 we have generated this set of rules by hand and checked both ourselves and 
by computer that TZ is correct. For K = 10 000, we have computed the values given in 
(TT|)- (fTU|) . see Figure HI and found out that 

SlMgHMgJ) = 0.299871 . 
K 

Lemma [3] now yields Theorem [IJ 

We have also constructed (with computer assistance) a correct set of 3607 rules such 
that Lemma [3] applied for K — 1 000 000 yields that every cubic graph has domination 
number at most 0.299309 n + 0(n/g) but we neither present nor claim this bound here. 
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+ + + ? -> 0000 o««««x 000000 

+ + X -> 000 o 000000 

+ + • -> 000 o • • • • o 000000 

+ + o -> 000 o • • • O + ? -> ©000000 

+ x ^ 0© o«««ox 00000© 

+ • ^ 00 o • • • o • -> ©00000 

+ o ^ 00 o • • • o o -> 000000 

X ^ o • • o + ? -> 000000 

• X + ? -> 0000 O..OX -> ©0000 

• XX -> ©00 o««o« + ? -> ©000000 

• X« -> ©00 o««o«X -> 000000 

• Xo -> 000 o««o«« -> 000000 

• • + ? -> 0000 o • • o • o -> 000000 

• • X -> 000 o • • o o + ? -> ©000000 

• •• -> 000 o««oox -> 00000© 

• • o -> 00© o • • o o • -> ©0000© 

• o + ? -> ©00© o • • o o o -> ©0000© 

• ox -> 00© o • o + ? -> 0000© 

• o« -> 00© o.ox -> 000© 

• o o -> 00© o • o • + ? -> ©0000© 

• + + ? -> 0000 o«o«x -> ©000© 

• +X -> 00© o«o«« + ? ©00000© 

• + • -> 000 o«o««X 00000© 

• + o -> ©0© o • o • • • -> 000000 
ox ^ ©0 o«o««o 000000 

o + + ? -> 000© o«o«o + ? ©00000© 

o+X -> 00© o.o.ox -> 00000© 

o + • -> 00© o • o • o • -> ©00000 

o + o -> ©0© o • o • o o -> ©0000© 

o o + ? -> ©00© o • o o + ? -> ©0000© 

o o x -> ©©© o • o o x -> ©000© 

o o • -> 00© o • o o • + ? -> ©00000© 

ooo -> 00© o.oo.x -> 00000© 

o • + ? -> ©00© o • o o • • -> ©0000© 

o • X -> 00© o • o o • o -> ©0000© 

o • • + ? -> 0000© o • o o o + ? -> ©00000© 

o««X -> 000© o.ooox -> 00000© 

o • • • + ? -> ©0000© o • o o o • -> ©0000© 

o • • • x -> ©000© o • o o o o -> 000000 



o 



Figure 3: The set of 79 rules used in the proof of Theorem CD 
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Figure 4: The numerical values of probabilities given by (|TT)-(|T0l) for the set of rules from 
Figure [3] and K = 10 000. 
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